Spin dynamics of dimerized Heisenberg chains 
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We study numerically the dimerized Heisenberg model 
with frustration appropriate for the quasi-lD spin-Peierls 
compound CuGeOs. We present evidence for a bound state 
in the dynamical structure factor for any finite dimerization 
8 and estimate the respective spectral weight. For the ho- 
mogeneous case {5 = 0) we show that the spin-wave veloc- 
ity Va is renormalized by the n.n.n. frustration term a as 
Vs = tt/2 J(1 - ba) with b ^ 1.12 
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Quantum ID spin systems may show a variety of in- 
stabilities. Of particular interest is the spin-Peierls (SP) 
phase due to residual magnetoelastic couplings [Q , which 
leads to the opening of a gap in the spin excitation spec- 
trum. The discovery of the spin-Peierls transition 
below Tsp — 14: K in an inorganic compound, CuGeOa, 
has attracted widespread attention. 

The spin-dynamics of CuGeOa is being studied inten- 
sively Above the spin-Peierls transition the one- 
magnon excitation spectrum, as measured by neutron 
scattering Q, forms a broad continuum. The form of 
the continuum is in good agreement with the continuum 
expected for one-dimensional antiferromagnetic Heisen- 
berg chains M; it is also seen in the magnetic Raman 
spectrum The physics behind the continuum in the 
one-magnon excitation spectrum lies on the fact that the 
elementary excitations of the one-dimensional Heisenberg 
chain are spinous (also called solitons or Bloch walls), 
each magnon being made up of two solitons. The result- 
ing "two-spinon" continuum has been seen experimen- 
tally also in other quasi- ID antiferromagnets [|o|. 

The exact form of the magnetic excitation spectrum in 
the dimerized state below Tsp is still being investigated. 
A two-spinon continuum is observed |^,^ with large spec- 
tral weight at the lower edge. The dispersion of the lower 
edge has been used to extract the spin- wave velocity [|j . 
Muthukumar et al. have pointed out recently that 
the experimental magnetic Raman intensity indicates a 
well defined magnon mode in the spin-Peierls state of 
CuGeOa and that the origin of this magnon could not 
be determined on the basis of the Raman scattering data 
obtained numerically for a single chain. 

The existence of a well defined magnon mode, i.e. a 
two-spinon bound state which splits of the continuum, 
has also been addressed recently by Uhrig and Schulz 



PH within an RPA approach. A recent neutron scatter- 
ing experiment p^ ] , has been interpreted as indicative of 
such a bound-state. 

In this context it is an important question to address 
whether other techniques can shed light on the magnetic 
excitation spectrum of the dimerized Heisenberg model. 
Here we present data for the dynamical structure fac- 
tor S{k,Lu) obtained by applying the recursion method 
|l3| to calculate the excitation spectrum of the hamilto- 
nian given below. This approximate method in particular 
gives very accurate results for the low-lying excitaions. 
We find evidence for a bound-state for any dimeriza- 
tion (5 > 0, independently of the amount of frustration a 
present. The finite-size corrections of the data is, for cer- 
tain k-values, small enough to determine accurately the 
position of the bound state and its spectral weight. We 
furthermore determine the renormalization of the spin- 
wave velocity in the homogeneous state as a function of 
the frustration. 

The magnetic properties of CuGeOs can be modeled 
by the ID spin-Hamiltonian 



H = jJ2[{'^ + Si-iy) S, • S,+i + aS, ■ S 
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where S is the dimerization parameter that vanishes 
above Tsp-, its zero temperature value has been estimated 
to be (5 = 0.03 The exchange integral J > and 

the amplitude for the frustrating n.n.n. coupling a have 
been estimated to be J « 150K ^ and a « 0.24 - 0.36 
p^ , p^ , respectively. 

The phase diagram of H in Eq. (|l]) has been calculated 
using the density-matrix renormalization-group method 
|p7| . For S — and a < etc « 0.2411, the ground state is 
gapless and renormalizes to the Heisenberg fixed point. 
The rest of the phase space is gapped. For 6 + 2a = 1 the 
exact ground state is known to be a valence-bond state. 

The dynamical structure factor S^k^uj) in the ground- 
state |0) of (U) is given by 

S{k,uj) - l(0|^"(-fc)l«, k)\^S{uj - {E^{k) - Eo)), 
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(2) 

where \n, k) and En{k) are the respective eigenstates and 
the eigenvalues of (|l|) in the subspace with total momen- 
tum k, |0) being the overall ground-state and Eq the 



1 



ground-state energy. The allowed values for the momen- 
tum k are k = 0, 27t/N, . . . , tt, for a finite chain of length 
N with periodic boundary conditions. S^{k) is given by 



S^k) 



1 



N 



J2 exp{ikx)S'^. 
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We have evaluated S{k,u!) for chains with N = 
12, 16, 20, 24 using an approximate scheme for the 
determination of the lowest-lying excitation energies 
En{k, N) — Eq{N) and the corresponding transition prob- 
abilities 1(7115^(^)10)1^ |13|. Using a recursion algorithm 
a set of orthogonal states is built starting with 5'^(A;)|0) 
[Q . Coefficients occuring in this procedure form a tridi- 
agonal matrix whose eigenvalues and eigenstates deter- 
mine the excitation energies and transition probabilities. 
(For further details and technichal limitations of this 
method see ref. Q). 

The dynamical structure factor is given, for any finite 
chain, by a sum over discrete poles. The weight of the 
individual poles, 

Wn{k) = |(0|5^(-fc)|n,fc)|2 

will go to zero in the thermodynamic limit, N oo, 
whenever the respective pole contributes to the contin- 
uum. The sum-rule 



Wnik) = Sik) 



(3) 



is here valid, where S{k) is the static structure factor, 

S{k) = {0\S^{-k)S%k)\0). (4) 

We are here interested in particular in the lowest pole, 
located at E{k) — i?o(fc) — Eq, where \n — 0,k) is the 
lowest state with energy Eo{k) in the subspace of total 
momentum k. We want to examine the question whether 
this pole is part of the continuum in the thermodynamic 
limit, or whether it evolves into a bound state, charac- 
terized by 



wo{k) 



const. 



oo. 



(5) 



In Fig. we present for a = the relative weight of the 
lowest pole contributing to S{k,u!), 



w'^^''\k) = wo{k)/S{k). 



(6) 



The data for Wq (fc) presented in Fig. (|l|) shows a mono- 
tonic decrease for (5 = as a function of chain length 
iV, indicating that the lowest pole will be part of the 
two-spinon continuum in the thermodynamic limit. This 
behaviour is what we expect for the homogeneous Heisen- 
berg chain. For the dimerized Heisenberg chain, 5 = 0.1 
in Fig. (|l|), we observe on the other hand for all fc < tt 
a monotonic increase of w^^^^ (k) with increasing chain 



length! This behaviour is also obtained when we intro- 
duce a finite amount of frustration a = 0.24 into the 
system, as can be seen in Fig. (^, where we compare the 
data for (5 = and 6 = 0.03, as apropiate for CuGeOs. 

We also observe in Fig. (|l|) and Fig. (||) that the finite- 
size dependence of the relative weight of the lowest pole 
is small enough in the dimerized state for k < 7r/2, that 
an estimate of the spectral weight in the thermodynamic 
limit is possible for k < Tr/2. Taken together we con- 
clude that the data presented in Fig. (|l|) and Fig. (||) 
clearly indicate for S > a finite value for wo{k) in the 
thermodynamic limit and therefore a bound state. 

Another interesting feature of the data presented in 
Fig. (|l|) and Fig. (||) is the large dependence of w^'^^^ (k) 
on the dimcrization parameter S. This observation is 
reminiscent to the large (5-dependence observed for the 
Raman spectral weight in a previous study |^. It is in- 
teresting to note that this behaviour is solely due to the 
(5-dependence of the matrix element wo{k) and that the 
static structure factor S{k), entering Eq. (||) is relatively 
size and ^-independent for most fc < tt, as can been seen 
form the data presented in Fig. (||). 

In Fig. we present for a = 0.24 the data for static 
structure factor S{k) and the energy dispersion E{k) = 
Eo{k) - Eo for both 5 = and (5 = 0.03 and for iV = 20 
(open symbols) and = 24 (filled symbols). We notice 
that the finite size effects are, in general, much smaller 
than those for the weight of the lowest pole presented 
in the previous figures and that the finite size effects for 
6 = 0.03 are smaller than those for (5 = 0. 

A relative large difference for ^ = and S = 0.03 oc- 
curs, for the data presented in Fig. for the "magnon 
dispersion" E{k), for k > 7r/2. This behaviour is a con- 
sequence of the coupling of the momenta fc and k + n 
by the dimerization. For any finite S > we have 
therefore the symmetry E{tt/2 + Afc) = E{7r/2 - Ak). 
In the homogeneous case, 6 — 0, one has instead that 
E{'k/2 + Afc) > E{Tr/2 - Afc) for any finite chain and 
tt/2 > Afc > 0. 

In Fig. (^) we present a study of the spin- wave velocity, 
Vs, as a function of frustration a and a range of system 
sizes A^ = 12; 16, 20, 24, and (5 = 0. Here we estimated 
Vs using the formula 



k=2-!v/N 



Strictly speaking we expect Vg to diverge like ~ AN/2tt in 
the thermodyamic limit for any a > ac « 0.2411 where a 
gap A opens in the spin-wave spectrum. This divergence 
will be very difficult to observe in finite-size studies as 
the gap is exponentially small for a larger but close to 
etc 
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The size-dependence of Vs is evident also for small val- 
ues of a. For a = the exact values for the spin- wave ve- 
locity is known to be 7r/2 J in the thermodynamic, given 
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by the intersect of the sohd Hne in Fig. with the y- 
axis. The finite size corrections do, on the other hand, 
actual increase shghty with increasing chain length N, at 
least for < 24 and a = 0. We have therefore decided to 
try for a linear fit of Vg as a function of a by demanding 
the fit to reproduce the Bethe-Ansatz value for a = 0. 
The corresponding fit (solid line in Fig. (0) is 



-J{l-ba), 



with b w 1.12 m 



In conclusion we have presented numerical evidence for 
the ocurrence of a discrete-pole contribution to the dy- 
namical structure factor for dimerized Heisenberg chains. 
This contribution to S{k,uj) may also be interpreted, at 
least in part, as a two-spinon bound state |]ll[| , it con- 
stitutes a well-defined magnon mode. It is interesting to 
note that such a bound-state would imply a peak in the 
two-magnon density of states and might therefore show 
up in the Raman-spectral weight. A previous numeri- 
cal numerical study has not found such a peak structure 
in the Raman spectral weight of dimerized Heisenberg 
chains. This discrepancy may be either due to strong 
magnon-magnon interaction effects in ID spin chains or 
due to matrix-element effects |^,|l^. 

We have furthermore presented data for the disper- 
sion of the magnon-mode and the static structure fac- 
tor, showing that the finite-size effects are small for these 
quantities. In addition we have estimated the renormal- 
ization of the spin-wave velocity for frustrated Heisenberg 
chains in the gapless region. 
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FIG. 1. The relative weight ?j;g'^°'' (fc) of the lowest pole 
contributing to the dynamical structure factor S{k, to) for 
a — 0, 6 = and S — 0.1 and for chains with 
N — 16, 20, 24 sites respectively. The relative weight ist given 
by ■w^''''\k) — wo{k)/ ^ Wn{k), where the w„{k) are the ab- 
solute weights of the poles contributing to S{k,ui). The lines 
are guides to the eye. 



ACKNOWLEDGMENTS 

This work was supported by the Deutsche Forschungs- 
gemeinschaft, the GraduiertenkoUeg "Festkorperspekt- 
roskopie" . 



FIG. 2. The relative weight w[^'^^\k) of the lowest pole 
contributing to the dynamical structure factor for a = 0.24, 
(5 = and 5 = Q.l and for chains with N — 16, 20, 24 sites 
respectively. The lines are guides to the eye. 
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FIG. 3. The energy of the lowest excited state with mo- 
mentum k, E(k), and the static structure factor, S(k), for 
a = 0.24 and both 6 = and S = 0.03. The open/filled sym- 
bols are the data for chains with A = 20 and A = 24 sites 
respectively. The lines are guides to the eye. 



FIG. 4. The spin- wave velocity, Vs, as a function 
of Q calculated for 5 = and chains with length 
N = 12, 16, 20, 24. Vs has been determined using the for- 
mula Vs = (E{2tv/N) - _B(0))/(27r/A), where E(k) is the 
ground-state energy for chains with total momentum k. The 
dahsed lines are guides to the eye, the full line is a fit. 
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